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LONG-TIME EXISTENCE OF CLASSICAL SOLUTIONS TO A 1-D SWELLING GEL 

M. CARME CALDERER t AND ROBIN MING CHEN* 

Abstract. In this paper we derived a model which describes the swelling dynamics of a gel and study the system in one- 
dimensional geometry with a free boundary. The governing equations are hyperbolic with a weakly dissipative source. Using 
a mass-Lagrangian formulation, the free-boundary is transformed into a fixed-boundary. Wc prove the existence of long time 
-solutions to the transformed fixed boundary problem. 



1. Introduction. In this paper, we study existence of long time C^-solutions of a free boundary problem 
modeling swelling of gels, in one space dimension. The gel is assumed to be a mathematical mixture of polymer 
and solvent. We assume that the gel is surrounded by solvent and that its boundary is fully permeable 
to solvent. The governing equations consist of the laws of balance of mass and linear momentum of the 
components and form a weakly dissipative hyperbolic system. The scaling leading to hyperbolic dynamics is 
motivated by the modeling of polysaccharide gels occurring in nature, for instance as motor devices of gliding 
bacteria (myxobacteria) . In general, hypebolic dynamics is a signature property of active gels encountered in 
models of the cytoskeleton of living cells and also in studies of biomimetic non-equilibrium gels [14] , [19] . 

We assume that a gel is a saturated, incompressible and immiscible mixture of polymer and solvent. 
Incompressibility in the context of mixtures refers to every component having constant mass density in the pure 
state. (In particular, it does not preclude the polymer from experiencing deformations with very large change 
in volume, as much as 500 percent, if ions are present). Immiscibility corresponds to the constitutive equations 
depending on the volume fraction, (f>i{= (j)) and 02, respectively, of the polymer and solvent component. The 
saturation assumption (also known in some literature as incompressibility) is the statement of the volume 
fractions adding to 1. The balance laws consist of the equations of balance of mass and linear momentum of 
each component together with the saturation constraint. The coupling between component equations takes 
place through the Flory-Huggins energy of mixing, the drag forces and the boundary conditions. 

For gels made of entangled polymers such as those used in biomedical devices, scaling arguments show 
that elastic ( [501 HU [53] ) and dissipative effects often dominate inertia, that enters the early dynamics only. 
Moreover, in [5] the authors assume that the dissipation of both the polymer and solvent is Newtonian. 

In this paper, we consider gels made of polysacharide networks, with dissipation parameters several orders 
of magnitude smaller (as much as 10~^) than those of the entangled polymer counterparts [22]. This motives 
us to keeping the inertia terms to study the intermediate time dynamics of the polysaccharide systems (as 
well as the early dynamics of the highly dissipative gels). 

The energy of the system is the sum of the nonlinearly elastic stored energy function of the polymer 
and the Flory-Huggins energy of mixing. We establish a dissipation law satisfied by solutions of the governing 
system in arbitrary space dimension. There are two types of boundary conditions complementing the elasticity 
and the Euler equation. These are of traction-displacement kind, for the gel equation (this being the sum of 
the elasticity and the Euler equation for the fluid) , and either scalar Dirichlet or traction boundary conditions 
for the the Euler equation of the fluid. These boundary conditions, also considered in earlier works on gels 
[55] , [5] , are motivated by the works of Doi and Yamaue [50] [5T] . They express properties of permeability 
of the gel interface to the surrounding fluid, and may range from impermeable (Dirichlet boundary) to fully 
permeable (traction boundary). However, these conditions have been found not be exact, from the point of 
view that they do not necessarily guarantee the balance of mass, linear momentum and energy at the interface 
between the gel and the external fluid. In fact, they differ from the balance laws at the interface by a term 
quadratic on the normal component of the difference of the polymer and solvent velocities, and proportional 
to the fluid inertia. Moreover, the assumption of viscous interface stress yields a fully dissipative energy law 

In this work, we consider a polymer separated from the surrounding fluid and confined to a strip domain. 
We assume that at time t — the polymer enters into contact with the fluid at the boundary points, causing 
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it to enter and swell the polymer. The requirement for the three dimensional governing system to have one- 
dimensional solutions determines the (pressure) Lagrange multiplier. Moreover, it follows from the equation 
of balance of mass, that the center of mass velocity is independent of space. This motivates us to rewrite the 
governing equations in terms of the relative fluid-polymer velocity, U = vi — V2. Another implication of the 
one dimensional geometry is that the volume fraction on the boundary is fully determined from the traction 
boundary condition, provided the fluid inertia on the boundary conditions is neglected. In this case, a problem 
originally formulated with traction boundary becomes a of Dirichlet type. 

Without loss of generality and to simplify the analysis, we assume that the initial domain (also used 
as reference configuration) to be the symmetric interval (—1,1). We use the mass-Lagrangian change of 
variables from gas dynamics, to transform the free boundary problem into one with fixed domain. We find 
that G'{(j)) -I- C/^ < is a necessary and sufficient condition for the governing equation to be hyperbolic, and 
so, it is necessary that G'{(j)) < hold. The function G((/)) corresponds to a one-dimensional stress, and the 
condition of negative derivative with respect to 4> expresses monotonicity of the stress with respect to the 
deformation gradient ^ where 0/ denotes the reference volume fraction of the polymer. Heuristically, 
if G' is bounded, the hyperbolicity condition is satisfied for sufficiently small relative velocity. Using the 
framework introduced in [2 in dealing with the hyperbolic initial-boundary-value problem, we obtain the 
local well-posedness of strong solutions for data satisfying the hyperbolicity and uniform Kreiss-Lopatinskii 
condition. We also consider the long-time existence of classical solutions. We use the method of characteristics 
to estimate the C^-norm of solutions. Due to the lack of strong dissipation of the system and the boundary 
damping, our estimate does not lead to the existence of global solutions. Instead, we prove that the existence 
time is at least of order O (| loge|), where e is the C^'^-deviation from the equilibrium. 

The rest of the paper is organized as follows. In Section[2]we set up the mathematical model that describes 
the swelling dynamics of a gel and derive the constitutive equations. We also formulate the one-dimensional 
system with boundary conditions. In Section [3] we transform the one-dimensional free-boundary problem into 
a fixed-boundary problem using the mass-Lagrangian change of variables, and prove the local well-posedness 
of strong solutions to the transformed problem. In Section [4] we prove the long-time well-posedness of classical 



2. General model of gel-swelling. We begin by setting up the model of the hydrogel postulated in 
21 E]- We assume that a gel is a saturated, incompressible and immiscible mixture of elastic solid and fluid. 
In the reference configuration, the polymer occupies a domain il G M''. The solid undergoes a deformation 
according to the smooth map 



We let = x(f2) denote the domain occupied by the gel at time t > 0, F C d^lt be the gel-solvent interface 
(in the case that T = dflt , the gel is fully surrounded by the solvent) , and let F = Vx denote the deformation 
gradient. We label the polymer and fluid components with indices 1 and 2, respectively. A point x S ilj is 
occupied by, both, solid and fluid at volume fractions 0i — (/)i(x,i) and 02 — 02(x, i), respectively. We let 
vi = vi(x,i) and V2 = V2(x, i) denote the corresponding velocity flelds. 

An immiscible mixture is such that the constitutive equations depend explicitly on the volume fractions 
4>i. We let Pi denote the mass densities of the ith component (per unit volume of gel). They are related to 
the intrinsic densities 7^ by pi = 7i0i, i = 1,2, which are constant in the case of an incompressible mixture. 

2.1. Governing equations. The assumption of saturation of the mixture, that is, that no species other 
than polymer and fluid are present, is expressed by the equation 



(In some terminologies, this is condition is also known as incompressibility). The governing equations consist 
of the balance of mass and linear momentum of each components as well as the chain-rule relating the time 
derivative to the gradient of deformation with the velocity gradient: 



solutions. 



(2.1) 



X = x(X,<), such that det(Vx) > 0, X e f7. 



(2.2) 



01 + 02 = 1- 



(2.3) 



dt 



^+V-(0,v,)=O, 



(2.4) 
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(2.5) 



Ik 



+ {vi-V)F^ (Vvi)i^, 



X e fit, i = 1, 2. Here fi represent drag forces between the components, and 71 is the Cauchy stress tensor of 
the ith component. Note that we have used the incompressibihty assumption 7^ = constant in the previous 
equations. The Lagrangian form of the equation of balance of mass of the polymer component is 



(2.6) 



where 



(f>i det F — 4)1 in 51, 

is the initial polymer volume fraction in the reference configuration. Note that if deti^ — 1, then 
1, so that no changes in the volume fraction of the gel correspond to no changes in volume fraction. 



Hence equation ( 2.6 1 is equivalent to equation (2.3 ) with i = 1. Adding equations (2.3 ), and using the constraint 



(2.2) gives 



(2.7) 



V • ((/)lVi + 02 V2) = 0. 



With this formulation, the final system of the governing equations are (2.4) and (2.5), subject to the 



constraints in (2.2), (2.6), and (2.7 



2.2. Boundary conditions. We now specify boundary conditions at the gel-water interface F. We 
observe that the equations of balance of mass are first order, scalar equations for 0i and (j)2 that require 
prescribing initial conditions only. Later, we will assume that the network part of the gel is elastic, with 
Newtonian dissipation. The latter will also be assumed for the in-gel and exterior fluids. This requires one 
vector boundary condition for each equation. There are two kinds of boundary conditions that we will discuss, 
in order to come up with a total of six scalar boundary conditions for the system: 

• statements of balance of mass, linear momentum and energy across F, and 

• relations expressing the permeability properties of the interface. 

Balance of mass across F. We denote v/ the velocity field of the fluid outside the gel network, which is 
assumed to be Newtonian. The equation of balance of mass of the fluid across the interface F is 



(2.8) 



(v/ - vi) 



)(v2 - vi) • n 



w, 



where n denotes the unit normal to F, pointing from the gel to the pure liquid. Note that the last relation 
defines the normal component, w, of the relative velocities. We also let ti and t2 denote a pair of orthonormal 
vectors perpendicular to n. Let v be a vector field on F. We use the notation || to denote the vector components 
tangent to F at a point. That is, we write 



= (v • n)n + Vj| = (v • n)n + ^ (v • ti)ti 



Continuity of the components of the velocity tangent to F. We assume the following no-slip condition 
(2.9) (vy - vi)|| = (V2 - vi)|| := q. 

Remark 2.1. The no-slip assumption equating the tangential components of the fluid velocity, relative to 
the polymer one, across F may admit a generalization of the form 



(2.10) 



(V/ - Vl)|| = 77„(V2 - Vl)|| 



where rjy > Q is a dimensionless parameter related to the surface viscosity. We will observe next that this 
brings tangential stress components in the balance of linear momentum equation. 

Balance of hnear Momentum across F. We next consider force balance across the interface F. Consider a point 
X g F at a time t, and let us observe this point in an inertial frame traveling with the same velocity as the 
polymer at point x and time t. 
travels at velocity wn 



Given equations (2.8) and (2.9) we observe that the water in the fluid region 
q and in the gel region travels with velocity -r^n + q. the change of mass of water 
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across the interface per unit time at point x at time t is given by 72|w|. As water comes out of or into the gel 
region, corresponding to the collapsing {w > 0) or swelling state {w < 0), respectively, there is the following 
amount of momentum change from the gel region to the water region per unit time: 



(2.11) 



72|w| [ wn - — n 



The force acting at the interface is given by: 

(2.12) {Ti+T2-T)n. 

The law of balance of linear momentum gives 

(2.13) 



(7i + 72 - T)n = 72|w|w ^1 - ^ j n. 
Note that this implies, in particular, that 

(2.14) (Ti + T2- T)n • t, = 0, ^ = 1, 2. 

Note that the difference in the stress across the boundary is normal to the interface. In particular, we have: 



(2.15) 

since q is tangential to the surface F. 



(rin + r2n-rn) •q = 



At this point, we have the boundary conditions (2.81, (2.9) and (2.13) on F, which together give us six 



boundary conditions. Mass balance and total force balance would provide the necessary number of boundary 
conditions if the interior of Qt were composed of a one-phase medium. Here, the interior of fit is a two-phase 
gel. We thus require three additional boundary conditions, assuming that all of the phases have bulk viscous 
stresses. This corresponds to specifying some condition that involves w and q. The appropriate forms for 
these boundary conditions will be discussed shortly. 



We now check that the boundary conditions (2.8), (2.9) and (2.13) lead to mass and momentum conser- 



vation. Let us first check that the amount of water is conserved: 



dt 



dx 



= J 02(v2 - vi) • ndS - y (vf - vi) • ndS 



=0, 



where we used (2.8) in the second equality. 



Now, we turn to momentum conservation. We have: 



d 
di 



/ (7101 Vi + 72(/>2V2) dx+ 72V/Cfx 

= j (7202 |(vi - V2) • n| V2 + Tin + 72n) dS 

-J (72 |(vi - V/) ■n\vf+ Tn) dS 

= J (72(V2 - v/)|w| + (Tin + T2n-T n)) dS 

= J (^72 - 1^ w\w\n + (Tin + T2n-T n)^ dS 

= 
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where we used ( 2.8 ) in the second and third equahties, ( 2.9 1 in the third equahty and ( 2. 13 1 in the last equahty. 



Balance of energy across T. Next, we consider energy conservation. The final form of our energy relation will 
lead us to possible forms for the additional boundary conditions we shall impose on our system. 

(2.16) 



d 
dt 



/ (vi(v-ri) + v2(v-r2))dx+ / vy(v-r)dx 

+ j Q7202||v2iP(vi - V2) • n - ^72||v/|p(vi - V/) • dS 

-/ ((Vvi)ri + (vv2)r2)dx- / {vvf)rd^ 

Jilt Jni 
+ ^ ((Tm) • vi + [T^n) ■ V2 - (Tn) ■Wf)dS 

+ j 07202||v2|p(vi - V2) • n - ^72||v/|p(vi - Vf) • dS. 



Let us evaluate the last two boundary integrals. Using (2.8) and (2.9) we have: 
(2.17) (rin)-vi + (r2n)-V2-(rn)-v/ 

:(TLn + 72n - Tn) • vi ' ' ' ^ 

1 



n - n • (Tn) w + (Tn - Tn) • q 



72^^ ( 1 - — I (vi • n) + ^n • ) - n • (Tn) ) w - (Tn) • q 



where we used (2.8) and (2.9) in the first equality and (2.15) in the second equality. On the other hand, 

1 



(2.18) 



272'/'2||v2|P(vi - V2) • n - ^72||v/|p(vi - Vf ) • n 



1 II 
=-272||vi 



™ 11? 1 II II? 

— n + q ^w + -72 vi +wn + q^w 

92 2 



1 



1 



W(l-:^)(vi-n)-U72 



1 



-72W w 



where we used (2.8|) in the first equality. We may go back to (|2.16|) to conclude that: 
(2.19) 



d 
dt 



(X (^'^^'^^"'^^"^ + ^7202|lv2||^^ '^^^ In': ^^2"'^^"^'^'') 
=- / ((Vvi)Ti + (Vv2)T2)dx- j {ywf)Td^ 

" /r ( (" • (^"^ - " • (I) " - (S) 

-^(Tn)-qdS. 

The last two boundary integrals denote the change in energy coming from the surface F. We would like these 
terms to be negative. One way to achieve this would be to let: 



w dS 



(2.20) 
(2.21) 



, 1 ■}\ I /T\ 1 I w 

T]±w^ [ n • (Tn) - -72W I - I n 



n 72 

\(P2 



V\\(l= (Tin) 11, 
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where ri± and ?7j| are positive constants and (7in)j| denotes the component of Tin that is tangential to the 
membrane. The above conditions provide the additional three boundary conditions we need on T. The 



boundary condition (2.201 depends quadratically on w and is physically reasonable only if w is sufficiently 
small. This difficulty will not arise if we neglect inertial terms and set 7^ = 0,i = 1,2. For most practical 



situations, inertial effects can be safely neglected. If we substitute (2.201 and (2.21) into (2.19), we have: 




c 2 



72||v/|pdx 



/ ((Vvi)ri + (Vv2)r2)dx- / (Vv;)rdx- f {r^^w^ + r^Hf) dS. 
Jilt -i^t "'r 

2.3. Energy dissipation and constitutive equations. The free energy density 4" of the gel consists 
of the elastic energy 'Wp{F) of the polymer and the Flory-Huggins energy of mixing yVFH(0i,02) [lO l 111 ) : 



(2.22) 
(2.23) 



Wfh(0i, 02) = a0i log 01 



KbT 

VmNi 



KbT 



62), with 

- 602 log 02 +C0102, 

KbT 



2V„ 



-X, 



where various parameters above are of the following physical interpretations: 

1. Kb is the Boltzmann constant, and T is the absolute temperature; 

2. Vm is the volume occupied by one monomer; 

3. A^i, N2 are the numbers of lattice sites occupied by the polymer and the solvent, respectively; 
4- X = x(0i: 02) is the Flory interaction parameter; 

In this way, the total energy of the gel is given by 



(2.24) 



nit 
■■=£p 



71 

T 



'1 vi 



72 



'2 V2 



* Ux 



72 I 



We now show that smooth solutions to the governing equations satisfy a dissipation inequality, and the 
dissipation inequality suggests the exact form of the Cauchy stress tensors for the ith component, %■ We 
define % as the sum of reversible stress T"^'^'' and viscous stress X*''''' for i = 1, 2: 



(2.25) 



(r) 



t: 



(v) 



Specifically, the dissipation inequality allows us to obtain the exact form of the reversible stress 7^ , the 

viscous stress T^^\ and the expressions for the friction forces fj. The derivation of the dissipation inequality 
uses a similar approach as in [5 , and we here follow the presentation from [5| . 



Theorem 2.2. (Dissipation Relation) Suppose that {vi,0i} are smooth solutions of equations [2.2)-(2.1) 
with boundary conditions (2.8), (2.9), (2.13), (2.20) and \2.21 ) on T. Let p denote the Lagrange multiplier 
corresponding to the constraint {2.1). Assume that the following constitutive equations for the stress tensors 
components hold: 



(2.26) 

(2.27) 
(2.28) 



T' 

'9. 



(r) 

J 

-{v) 



-02PI, 



902 



7-^"^ =77,Z?(v,) + A*«(V-v,;)I 



with i = 1,2, where rji > and fii > denote the shear and bulk viscosities of the ith component, respectively, 
and D(v) := ^(Vv + Vv"'") is the symmetric part of the velocity gradient. Suppose that the solvent-polymer 
friction forces are given by 



(2.29) 



fl = pV0i - k(vi - V2), f2 = -fl, 
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with K > denoting the gel permeability coefficient. Then the following dissipation relation holds: 



d£ 

Itt 



(2.30) 



■i 



where n represents the unit outer normal to fit. 

Proof. Without loss of generality, assume that 71 = 72 = 1- We observe that the boundary of the gel is 
determined by that of the polymer network and moves with velocity vi . We first calculate the time derivative 
of £. Using the Reynolds Transport Theorem, and taking into account that the boundary d^t moves with the 
network speed Vi, we have 



d£p 



V- ( vi^llvil 



9* 
~dt 



V • (vi0iWp) + V • (viWph) 



Next, expand ^ in terms of 0i, 7 ^ using the chain rule to obtain 



V-(v,f||v2|P 



d£p 



'1 II l|2 , II ||2 

9^2-""^" +yl|v2|| 



V - ( vi^llvil 



fit 



Wp 



9Wp 



5l,t 



[>2.t 



, dWp 



Ft 



(V • vi)0iWp + vi • V(0iWp) + V • (viWph) 



dx. 



dx 



Equation (2.3 1 of balance of mass, = —V • (0iVi), for i = 1, 2 yields 



d£p 
"dT 



1 II l|2 I ^2 II ||2 



V-(vi^||vi|p 



V-(v,f||v2|P 



- Wp + 



901 



V • (</)iVi) 



d(f>2 



V • (02 V2) + 01 



awp 



dx 
dx 



(V • vi)0iyVp + vi • V(0i>Vp) + V • (viWph) 



dx. 



Gathering the terms involving Wfh a-nd using the incompressibility of mixture constraint in equation (2.2) 
gives the following: 



V • (01 vi) 



awp 



-V • (V202) + Vi • VWpH + WphV • Vi 



V- (^0ivi)( 



aWpH ^awp^ 

901 902 

sWpH aWpH 



V a0i d02 



902 901 



Wp 



(V-vi 



WphV • vi 



Substituting it into the expression for gives 



d£p 
~dt~ 



'1 II l|2 , 'P211 ||2 
„ 1 — Vi H V2 

5t V 2 " " 2 " " 



V-(vi^||vi|p 



9WpH 9Wp 



WpH 



V-(vif||v2|P 



(V • vi) + 0ivi • VWp 



I dx. 
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Next, we apply the chain-rule relation between material and spatial time derivatives, 

DF 

to obtain the following identity: 



(2.31) 



Ft + (vi • V)F = (Vvi)F, 



dF 

^ dF 
dWp 

^ dF 



bl^^l ■ VWp 



: Ft + (/)ivi 



dWp 
dF 



: Vf] 



dWp 

: {Ft + VI • VF) = ct>i : VviF. 



The expression for ^^jp simplifies to 



d£p 
~dt~ 



^ yllvill +y||v2|| 



+ 



: ((Vvi)F) + 



dct>i 



+ >Vp 



dx 



(V • vi) } dx 



In the final steps, we address the kinetic energy terms. Using the equations of balance of mass and linear 
momentum of for i = 1, we calculate: 



d^ 
dt\2 



iiiviin+v- vi-<^iiivi 



= ^lVi •Vi,t+ ^||vi| 



dt 



V(||vi| 



<^ivi • vi,t - ^llvill^V • {(t)ivi) + ^llvill^V • (01 vi) + .^ivi • (Vvi)^vi 



~15t 



= vi-(v-ri + fi) 



Here, Mi denotes the material time derivative. Likewise, 



Dt 



d (\ 



dt\2 



■2||v2f ) +V- ( Vil(j)2\W2 



= (/>2V2 • V2,t + -||V2| 



dt 



V-(</>2Vl)) +^<^2Vl-V(||v2f) 



= (j!)2V2 • V2,t - ^||V2||^V • (?!>2V2) + ^||V2||^V • (h^l) + h^l ■ (VV2)'^V2 
= (f)2V2 ■ V2,t + (l>2^2 ■ (VV2)^V2 + (^2^1 ■ (VV2)^V2 
-^f>2V2 • (VV2)'^V2 + ^||V2||^V • (9i2(vi - V2)) 

= '/>2V2 • + (/'2(VV2)^V2 • (Vi - V2) + ^||v2|pV • ((/>2(vi - V2)). 



Integrating the last term by parts gives 



/ 



|(i*,llv,ll" 



V- ( Vi-02||v2f 



dx 



/ 



<^2V2 • + ?!'2(VV2)'^V2 • (Vi - V2) + ^||V2||^V • ((A2(vi - V2)) 



dx 



= / 02 V2 • dx 



2'?!'2||v2|P(vi - V2) • ndS. 



2-8 



Plugging in the Lagrange multiplier we have that, for i — 1,2, 

-- I y"(v» • (V • 71) + f. • V, +pV • (0,v,)) dx 



dEp 



0(1)2 d(j)i 



(V • vi)] \ dx 



+ ^ / <?!)2||v2|p(vi - V2) -ndS. 
Integrating by parts on the term V • 7i for i = 1, 2 gives, 

/ V'l-Vvi • Ti - p0i(V • Vi) + fi • Vi - pV(/)i • Vi) dx 



d£p 



(2.32) 



^i^:(Vvi)F- 



TiVi + -'?ii2||v2|P(vi - V2) 



•ndS. 



(V • vi) ^ dx 



The time derivative is treated in a similar way giving 



(2.33) 



dSs^ 
dt 



Vvy • r dx 



'^v/ + ^^(vi-v/) ) -nd^. 



Therefore the conclusion follows by adding up (2.32) and (2.33), and using Ti as in (2.26) and (2.28), fi 
as given by (2.29), and all the boundary conditions. 
□ 

We further assume that the polymer is an isotropic elastic material, that is, 

(2.34) Wp{F) = Wp(/i,/2,/3), h = trC, [tr2(C) - (tr Cf] , h = det C, 
where C — F^F. Here we take Wp{F) to be of the following form 

(2.35) Wp(F) = (11 -c) + ao{i;^ - 1) + /3o/| + 

where c, ao, Po, /3i > 0, r > 1 and s, q > I. 

As for the Flory-Huggins mixture energy, we adopt the interaction equation from Horkay et al. |13j as 
follows. 

(2.36) x{(f>i,(f>2) ^Xo + Xl(|>l+X24'l■ 
Fm:t]lCT computation gives the following results for the reversible stress tensors 

r/''^ = 2sr-^FF'^ + [-aor(det Fy + /3o det F + /3ig(det Ff] I 

-7pr^Xa<l>2 + -TTTy- log <pl + + 2xi(j)l(l)2 + 3X2(f>l<p2 



(2.37) 



(KbT ^KbT, , ^ KbT 

^-7^X001 + -TTTr- log ^2 + TTTT- 



V 2v; 

KbT 



N2Vm 

KbT 



N2Vm 

. KbT 



+ Xi<l>l + X2(f>l 



7,7^XoM2 + 01 log 01 + h log 02 



+X10?02 +X20?02 



(2.38) 



r 
'2 



(r) 
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2.4. A new field of unknowns. Now wc will formulate the effective governing equations in terms of 
the center of mass velocity V = (j>iVi + 4>2^2 and the diffusion velocity U = Vi — V2. Hence the new field of 
unknowns is 



{V,U,0i,p,i^}. 



The total stress of the system is 



(r) 



T 



(r) 



6i02U(g) U, 



(2.39) 
(2.40) 

(2.41) 

(2.42) 
(2.43) 



In this way, in Eulerian coordinates, the governing equations (2.2)-(2.7) become 
d(j> 



dt 
dV 

Ik 
dlJ 

Ik 



i + [(V + (1 - 0i)U) • V] 01 + 0iV • [V + (1 - 0i)U] 

(V • v)v = V ■ r, 



0, 



(1 - 20i)(VU)U - (U «) U)V0i + (VV)U + (VU)V 



1 



-V-Ti 



1 



/3 



-U- 



AV01 



<P1 1-01 " 0l(l-0l) ' 0l(l-0l)' 

Ft + [V + (1 - 0i)U] • VF = V [V + (1 - 0i)U] F, 

v-v = o. 



The first and second equations give the balance of mass and linear momentum for the mixture. The third 
equation can be interpreted as giving the evolution of the microstructure of the gel. Equation (2.42) is a 
version of the chain rule relating time derivatives of F with velocity gradients. This equation is required in 
mixed solid-fluid systems |18J . The last equation looks like an incompressibility condition of the mixture and 
is from the balance of mass and the constraint 0i + 02 = 1- 

In application to modeling gliding behavior of bacteria by polysaccharide swelling, the problem can be 
thought of as in one dimensional fl2\. We consider the gel occupying a strip domain {(cc, y, z) : —L < x < L} 
for some L > 0. Now the fields of problem become 

(2.44) V= («(x,i),0,0), U = (u(x,i),0,0), 0i = 0(a;,t), p = p(x, y, z, t). 
The deformation gradient is 

(2.45) F = diag(detF(a;,i),l,l). 



From (2.61, (p {x, t) det F{x,t) = 0/. We can also write down the stress tensors (2.37) and (2.38) in the 
one-dimension formulation 



2s I 



b^j(j)'^+2Y 'diag I 



.1,1 



KbT KbT 

Xo(l - 0) + ^TT-f^ log 



, , KbT , 



(2.46) 



2Vr 

( KbT 
KbT 



Xo0- 



KbT 



log(l 



(2.47) 



21^™ 
+Xi0'(l 

T2= -v{i - 4>)l- 



Xo0(l 



N2Vm 



)l0g( 



KbT 

N2Vm 

KbT 

' N2Vm 



+ 2x10(1-0) 

+ X10^+X20^ 

(1 - 0) log(l - 



/3ig0?0-9I 
i) -I- 3X20' (1-0) 



X20'(l 



The second and third equations of (2.40) indicate that p = p{x,t) is independent of y and z. Moreover, 
equation (2.43) together with the first component of equation (2.44) gives v = v{t). Prescribing the initial 
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condition i;(0) = leads to v{t) = for all t > 0, provided that V • T = holds. The latter determines p in 
terms of </> and u, up to a constant. In this way we arrive at the following system for (/) and u 



(2.48) 



where 



(2.49) 



0t + [0(1 - 0)u], = 0, 
'1 



-«^(l-2<^)-G(0) 



/3u 



G(0) 



j^i;riog(l 
N2V„, 

^ ^2 



j^Briog((/)) 



1 Vm 



(f)] + 202 



-(g-l)/3i0?r''-(l + r)ao07>'- 



2X10 + 3(xi-X2)0'+ 4X20' 



2.5. Boundary conditions. We assume initially, the polymer occupies the domain f^o — {{x,y,z) : 
—L < X < L}, and the solvent is in |x| > L. At a later time t > 0, the gel occupies the region fit — {{x, y, z) : 
—Si{t) < X < S2{t)}, where x = 5'i,2(i) are the positions of the interface between the gel and the pure solvent. 
Therefore inside the polymer region —Si{t) < x < S2{t) equations (2.48) hold for and u. 

The general boundary conditions are described in Section |2.2[ Here we further simply the problem by 
ignoring the inertial effects from the system and hence the stress balance (2.131 now becomes 

(2.50) {Ti+T2)n = Tn, on T. 

The other boundary conditions we impose here for our problem characterize the degree of permeability of 
the interface (see also, for instance, [HI [Midi])- Again with no inertial effects, (2.20) becomes 

'72^ 



(2.51) 



rjj^w — n • (Tn) — n ■ 



n. 



We now classify the boundary permeability in the following way 
(1) The interface is fully permeable if ?7_l =0. Thus 



(2.52) 



n • (Tn) — n ■ 



n = 0, on r. 



(2) The interface is impermeable if r/_L = oo. In this case w — that is 
(2.53) vi • n = V2 • n, on F. 

(3) The interface is semipermeable for ?7_|_ g (0, oo). 

In this paper we will consider the fully permeable interface. Hence from (2.501 and (2.52) we know that 

02 



In one dimension, plugging in (2.461, (2.47) and the above we obtain 







2s 



(2.54) 



KbT 
_ ( KbT 
{KbT 



0f0" 



Xo(l-0) + 



KbT 



Xo0- 



log(l 



log( 





KbT 



V 2F„ 
+X10'(1 



Xo0(l 



KbT 

N2V,n 
X20'(l-0)}. 



) log 4 

on r. 



KbT 
N2V^ 
KbT 



+ 2x10(1 -0) + 3x20'(l - 

+ X10^ +X2</ 

(l-0)log(l 
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Therefore we can determine the saturation value (p ~ (p* at the interface. In other words, the locations of the 
interface x — —Si{t) and x = S2{t) are determined by the saturation value (f)* , i.e. (j){t, —Si{t)) — 4>{t, S2{t)) — 

The kinematic boundary condition asserts that the interface moves with the speed of the polymer, which 
means -S[{t) ^ [I - 0(t, -5i(t))]u(t, -S'i(t)) and S'^(t) = (1 - (f>{t, S2{t)))u{S2{t),t). Therefore we have 
obtained the initial and boundary conditions as follows 



(j>{x,t) ^ (j)*, a.t X = -Si{t),S2{t) 
Si{0)^L, S[{t) = -[I ^ ^it,~Simu{t,-Si{t)) 
S2{0) = L, = [1 - ^{t, S2{t))Mt, S2{t)) 

(j){x, 0) ^ (j)", u{x, 0) = for - L < X < L. 



(2.55) 

3. The transformed problem. In this section we aim to setup a fixed-boundary problem associated to 
(2.48) and (2.55) and establish the local-wellposedness of strong solutions. 

To transform the free boundary to a fixed boundary, we perform the following change of coordinates 



(3.1) 



y 



Slit) 



(z, t) dz, T = t. 



S it'] 

Because dz gives the total mass of the polymer, we may normalize that to be 1. In this way the free 

domain {—Si{t), S2{t)) becomes the fixed domain (0, 1). Therefore the free boundary problem now turns into 



(3.2) 



+ Cjy^il - (j))Uy - (f)'^(f>yU = 0, 



_ -/3m 

y ~ </.(!-<<.; 



<^{y,T) = 0*, for y = 0,1 

0(2/, 0) = m(2/,0)=u", forO<2/<l. 



To further rewrite the system, we let ■;/; = 1 /0 and /(s) satisfy that /'(s) = sG'{s). Then let F{s) — f{l/s) 
In this way the above system becomes the following initial-boundary value problem 



(3.3) 



/ 



- ( 1- 

2 



U 

V 



1 





-/3uijj 



(1- 



(0,1) X (0,T) 



ip = V'*, at y = 0, 1, 



v. 



for < y < 1, 



where ip* = l/ip* and -0*^ = 1/0°- The gradient matrix is 

A{ip,u) = 



with eigenvalues 



(3.4) 



V 



02 



V;3 



u 



\l,2{i^,u) 



u^+G'{l/ij) (1-V) 



and the corresponding left and right eigenvectors are 
(3.5) L 



iMM={^T^^ 1-^ '^J 
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and 



(3.6) 



/ l-i; 



V 1 / 

In the range of physical parameters corresponding to semi-dry polymer we have < </) < 1, hence 1 < ij). 



Thus the system (3.31 is hyperbolic if 
(3.7) 



Hence G'{(j)) < will be needed to guarantee hyperbolic! ty of the governing system, and therefore is a 
requirement for the propagation of the swelling front towards the solvent region. It turns out that this 
condition is satisfied for polymer data (see Fig.l). However in the case of polysaccharide data, there may 
be multiple quantities (pc such that G'{(t>c) = (see Fig. 2). This may be interpreted in terms of the onset 
of deswelling, observed in bacteria motility phenomenon [12j : it may also be associated with volume phase 
transitions observed in systems with a small elastic shear modulus [IB] . 





Fig. 3.1. G and G' for polymer data 




0.4 0.5 0.6 

one — dimensional stress G{(p) 




Fig. 3.2. G and G' for polysaccharide data 



parameter 


iVi 


N2 


q 


s 


r 


a 


/3i 


0/ 


Xo 


Xi 


X2 


polymer 


1000 


1 




6 


1.25 


0.001 


20 


0.05 


0.467 


0.593 


-0.42 


polysaccharide 


1000 


1 


2 


0.6 


1.25 


0.001 


0.002 


0.05 


0.446 


0.106 


-0.02 



parameter values 



my- 



In concern with the local- wellposedness of system (3.3), we need to check the following conditions (see 



CI. Non-characteristic condition. The matrix A^ijj, u) is non-singular for (-(/;, u) in a certain proper domain 
M. 
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C2. Normality. The boundary matrix B is of constant, maximal rank and 



R^='kevB ® E'{A{i}:,u)) ai y = I 
= kerB©£;"(A(-0,it)) at y = 0, 

where E^{A{'ip,u)) is the stable subspace of A{'4>,u) and is the unstable subspace of 

Note that in our problem, the boundary matrix is 



B 



1 




C3. Uniform Kreiss-Lopatinskii (UKL) condition, for all {ip, u) £ M there exists C > so that 
(3.8) IIV^II < C||By|| 

for all V in the unstable subspace of (■(/;, u) at y = and for all V in the stable subspace of 
A~'^{ip,u) at y = 1. 



Consider now our problem (3.3 1, the first two conditions (CI) and (C2) are equivalent to 
(3.9) Ai(i/;,'a) < < A2(V',u), 



which is, from (3.4 1 
(3.10) 



< LJ^G'(1/V). 



Notice that since > 0, (3.10) implies that G'{\/^) < and moreover the hyperbolicity condition (3.7 1 



As for the third condition (C3), the stable and unstable subspaces of A ^(i/j,u) is spanned by i?i and R2 



(as defined in (3.6|) respectively. Hence the UKL condition (3.81 is satisfied when there exists a 7 > such 
that 



(3.11) 



■01 



1 - -0 



7/2 +G'(1/V') 



>7- 



In order to establish the wellposedness of the initial-boundary value problem (3.3) in some strong Sobolev 



space i7™ with m > 2 being some integer, the data should satisfy the compatibility condition 



for all p e {0, 1, 
(3.12) 



dPr{y,0)=df^{y,0), at y = 0,l, 
— 1}. Since "tp* is some constant, the above condition is simply 

af V'(y, 0) = 0, at 0, 1, for all p G {1, . . . , m - 1}. 



We can now state our local- wellposedness result. 

Th eorem 3.1. fE/Ifm > 2 is an integer, t hen fo r all (V^",!*") S iJ™+i/2([0, 1]) x 7J"+i/2([o, 1]) satisfying 
(3.1(f^ , (3.11), and the compatibility condition (3.12), there exists T > such that the problem (3.3) admits a 
unique solution u € i/'"([0, 1] x [0,T]). 



4. Long time existence of classical solutions. The Cauchy problem of system (2.48) was discussed 
in [3]. The authors showed that the Cauchy problem is L^-stable, and the source term is merely weakly 
dissipative. However, due to the result in j7|, the global existence of BV solutio ns c an still be obtained. 

In concern of classical solutions to the initial-boundary value problem of (3.3 ), the short time existence 



can be established using the general approach introduced in Chapter 4 of [IT. To obtain large time of global 
solutions, one needs to control the C^-norm of solutions. It turns out that if the systems exhibits strong 
enough dissipation or boundary damping then solutions can be extended globally in time (see, for instance 
[15j). However in our case, the above two types of damping are both weak. What we obtain is the large-time 
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existence and uniqueness of solutions to the system (3.3), provided that the data are close enough the 
equilibrium state {ip*,0). 

Let r] = t/j — ijj* . To simplify the notation we still choose {x,t) as the space-time variable. Then (3.3) 
becomes 



(4.1) 



B 



+ P{t],u) = 0, in (0,1) X (0,T) 



at a; = 0, 1, 



t=o 



where 



A{r],u) 



- if;* 



for < x < 1, 



(r/ + t/i* 



V 



+ G'{l/{r^ + r)) 







(?? + V^*) - 1 



Now we state our main result. 
Theorem 4.1. Suppose that tp* satisfies 



(4.2) 



G\i/r) > 



and the -compatibility conditions of the initial and boundary data 

' ??0(0) =??0(1) ^0, 

(4.3) \ 

I (^*)2^" ^ 
hold. Then for any Tq > 0, there exists an s > so that if 

||(»y",w")||ci <£, V < x < 1, 



then system (4--1) admits a unique classical solution for t G [0,ro)- 

Proof. As is pointed at the beginning of this section, the local-in-time wellposedness of solutions can 
be proved using the idea from [T7] . For the time being, suppose that on the existence domain of solution 
ri{x,t),u{x,t) we have 



(4.4) 



\{-q,u){x,t)\ < eo, 



where eo > is a suitably small number so that ( |3.10[ ) is satisfied for I77I, \u\ < Eq. Because of condition (4.2) 
and continuity, such an Sq exists. 

To get the large-time existence of solutions, it suffices to prove that we can choose Sq > small enough 
so that for any fixed e with < e < Eq, there exists some 6 = S{e) > small such that if 

(4.5) 11(77°, u°)||ci <S, VO<x< 1 
holds, then on the whole existence domain of the solution (rj, u) we have 

(4.6) \\{r],u){-,t)\\ci <s, yt>0. 

First we perform a diagonalization. Let 

'Ji = L^{r],u){r],u)^ {i = 1,2), 
jji^ Li{ri,u){dxri,dxuY (i = l,2), 
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where Li is the i-th left eigenvector defined in (3.51 with ij; replaced by 77 + ijj* . It is not hard to see that 
V — {wi,V2) and w — (ioi,z«2) satisfy the following system of diagonal form (cf. Chapter 3 in |15j ) 

2 2 
dtVi + XidxVi + k{vi + U2) = ^ CijkVjVk + ^ dijkVjWk 



(4.7) 



for i — 1,2, where 



2 2 

j,k=i j,k=i 



2(7/^* - 1) 



CijkT dijk, Cikj and dy-fc are continuous functions of {ri,u), Xi(t,x) = Xi{ri,u) (cf. (3.4 1 with ip replaced by 
rj + ijj*). Denote the quadratic parts by 

2 2 

Qi{v,w)= ^ CijkVjVk + ^ dijkVjWk (i = l,2), 

2 2 
Qi(u, w)= ^ CijkWjVk+ ^ dijkVjWk (i = l,2). 

The initial condition for w and w is obviously 

v,\t=o = v^ = L,{v',u°)iij°,u'^)^ [i = 1,2), 
u;,|,=o = = L.(7?o,wO)(a.77°,5,uO)^ (z = 1,2). 



Boundary condition for v: 



At X = 0, 1 : B 



L2 



Thus 

(4.8) At X = 0, 1 : vi= V2. 

Differentiating with respect to t, we get the boundary condition for w: 



= B 



= B 



-BA{ri u) 



L2 



-P{v,u)\ ^-BA{rj,u) 
-1 

Wl 
W2 



Therefore 
(4.9) 



At a; = 0, 1 : wi = W2. 



Therefore in order to prove (4.6), we only need to show that the same estimate holds for ||(i;,k;)||(7o: 

(4.10) \iv,w){x,t)\<e. 
provided that the initial data is small: 

(4.11) |(t'°,w°)| <S, V < X < 1. 
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From (4.4 1, we may also assume that on the whole existence domain of the C° solution {v,w), 
(4.12) \{v,w){x,t)\<SQ. 
Let 

Am 



A. 



Ti = 1/K 



min{|A,(t,a;)| : < t < T, < a; < 1, i = 1, 2}, 
: max{|Aj(t,a;)| : < < < T, < x < 1, « = 1, 2}, 



T2 = 1/A„ 



Since Eq is chosen so that the hyperbolicity condition (3.10) holds, we see that A,„in < and A,„ax > and 
hence Ti and T2 are well-defined. 



By continuity we know that we can certainly pick some S > small such that (4.10) holds on some time 



domain. Hence to prove (4.10), it is only necessary to show that there exists some Eq > so small that for 



any fixed T > 0, if ( |4.10[ ) holds on the domain D{T) = {(t, x) : < t < T,0 < x < 1}, then it still holds on 
D{T + Ti), provided that the C° solution {v, w) exists on such domain. 
For this purpose, let 



(4.13) 



Vi(t) = max \vi(x,t)\, Wi(t) — max \wi(x,t)\, i = l,2, 

0<2;<1 0<2;<1 

Ui{t) = max{Vi{t),V2{t)}, U2{t) = ma.x{Wi{t),W2{t)}. 



Suppose that the solution {v, w) exists on D{T + Ti) and let ^ = /^(t; x, t) be the i-th characteristic 
passing through a point (x, t) G D{T + f) with T <t <T + Ti. Then 



(4.14) 



dr 



fi{T\X,t) = \i{T,f,{T;X,t)), 

^t: fi{t;x,t)=x. 



From (3.9) we know that for vi{x,t) there are two possibilities: 

(1) The first characteristic ^ — /i(t; x, t) intersects the interval [0, 1] on x-axis with the intersection point 
(/i(0;a;,i),0), see Fig. 3. 



t 



Fig. 4.1. 



X 



Integrating the first equation in (4.7) along the first characteristic ^ = /i(r;x,i) we get 



vi{x,t) 



v°{MO,x,t))- / Ke-<'-^\2{MT;x,t),T)d7 



(4.15) 



e-^^'--'>Qi{v,w)ifi{T;x,t),T)dT. 



4-17 



From (4.11) and since (4.101 holds on D{T), we have 

/■* /■* ^ 

e">i(x,t)|< (5 + (e^^ - l)e + Ce^ + / Ke^'VaMdr + C / ^ K^^^t) + W^^^^)^^ 

< (5+(e'='^-l)£ + Ce2+ / Ke^Wi{T)dT + C f e'''' {Ufir) + Ui{T))dT, 



(4.16) 



where , and Ui were defined in ( 4.13| ) . 

(2) The first characteristic ^ = /i(r; a;, <) intersects the boundary x = 1 at point (1, ti{x, t)) where {ti{x, t)) 
satisfies 



Obviously we have 
(4.17) 



fi{Ti{x,t);x,t) = 1. 



T2 >t~TAx,t) > 0. 



In case (2), for the second characteristic ^ = f2{T; l,Ti{x,t)) passing through {l,Ti{x,t)), there are still 
two possibilities: 

(2a) This second characteristic intersects the interval [0, 1] on the x-axis with the intersection point 
(/2(0;l,ri(a;,i)),0), see Fig. 4a. 




(l,ri(x,t)) 









(l,Ti(a;,t)) 


(0,Ti2(x,f)) 





X (a) 
Fig. 4.2. 



X (b) 



(2b) This second characteristic intersects the boundary a; = at the intersection point (0, ri2(x,f)) (see 
Fig. 4b), where 

Ti2{x,t) = r2(Ti(x,t), 1), 

in which, T2(x, t) stands for the i-coordinate of the intersection point of the second characteristic ^ — /2(t; x, t) 
passing through a point (x, t) with the boundary x = 0: 



We have 
(4.18) 

Therefore 



Noting (4.17) we get 



h{T2{x,i)]x,t) = 0. 
T2 > i-T2(x,<) > 0. 

T2 > Ti{x,t) - Ti2ix,t) > 0. 
2T2 >t-Ti2{x,t) > Ti. 
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Therefore for T < t < T + Ti , 

(4.19) 0<Ti2(x,t) <T, T-Ti2{x,t)<2T2. 



In case (2a), using boundary condition (|4.8 ) we get 



e«*i;i(a;,t)= e«^i«i(l, n) - f ne''^ V2{h{T;x,t),T)dT - ( e^^Qiiv,w)ifi{T;x,t),T)dT 

- e'^^ii'2(l,ri)- / Ke^^V2{fiiT;x,t),T)dT- f e''^ Qi(v,w){h{T-x,t),T)dT. 
J n J Tx 



Integrating the second equation in (4.7) along the second characteristic we have 



e"^^«2(l,Ti) = z;°(/2(0;l,Ti))- / Ke^^vi{f2{T;l,n),T)dT 

Jo 

- f \''^Q2iv,w)if2iT;l,n),T)dT. 

Jo 

Combing the above equahties we get 

e^'vi{x,t)^V°2{f2{0;l,Tl))- f \e^^Vl{f2{T;l,Tl),T)dT- [ Ke^^V2{fl{T;X,t),T)dT 



JO Jti 

\^^Q2iv,w){f2iT;l,Tl),T)dT~ f e'^^ Ql{v , w){h{T- X ,t) , T)dT. 
Jti 



We have two cases: 

(i) Ti < T. Then it is easy to see that we have the same estimate as (4.16) 

(ii) T > T. Then 



e'''\vi{x,t)\< 5 + (e'^^ - l)e + Ce^ + J \e''^Vi{T)dT + Ke'^^T/2(r)dT 

2 



+C / e--Y.^,Hr) + Wi{T)dT 



(4.20) < S+ie"^ -l)e + Ce^ + f Ke^Wi{T)dT + C f e'^^(C/f (r) + C/|(r))dT, 

Jt Jt 



which is again, the same as (4.16) 



In case (2b), using the similar idea we obtain 



e'''vi{x,t)= e«^iwi(l,ri)- / Ke^^V2{fi{T;x,t),T)dT ~ / e^^Qiiv,w)ifi{T;x,t),T)dT 



e^^^V2{l,n)- f\ 

Jti 

e""^^W2(0,ri2)- / \e^^vi{f2iT;l,n),T)dT- f 

Jti-? J t^ 



Ke'^''v2{h{T\x,t),T)dT - / e'^^Qi{v,w){fi{T;x,t),T)d7 



Ke'^''v2ifi{T;x,t),T)dT 



Q2{v,w){f2{T;l,Ti),T)dT - / e""^ Q i{v , w){fi{T; X , t) , T)dT. 



From (4.19 ) we know that < T12 < T. Hence we can again spHt case (2b) in the foUowing two possibihties: 



4-19 



(i) n < T. Then the estimate for is 



\vi{x,t)\< 5 + ie"'^ - l)e + Cs'' + j ne''^V2{T)dT + C f e^^^V^ir) + W^{T)d7 

(4.21) < S + {e^^ -l)s + Ce'^ + [ Ke^^Ui{T)dT + C j e"-^ {U^{t) + U^{T))dT. 

Jt Jt 

(ii) ri > T. Then 

e'^>i(a;,i)|< 6 + {e"'^ - l)e + Ce^ + J \e''''Vi{T)dT + J Ke'^^FsM^r 
"'^ i=i 

(4.22) < 6 + {e"^ - l)s + Ce'^ + [ Ke''''Ui{T)dT + C [ e'^^(C/f (r) + Z7|(T))dT. 
Hence we always have 



(4.23) e''>i(a;,i)| < <5 + (e''' - l)e + Ce^+ / Ke''''Ui{T)dT + C e""^ {U({t) + Ui{T))dT. 

JT JT 

In the same way, we obtain similar estimates for \v2\^ \ wi\ and |w2|- Therefore we have 



(4.24) 

Now let 
(4.25) 



e''*f/i(x, t)< 5 + {e"'^ - 1)£ + Cff2 + Ke''Wx{T)dT + C e'^^Ufir) + Ui{T))dT, 
e''*U2{x, t)< S+ (e^"^ - l)e + Ce^ + Ke«^?72(T)dr + C e«^(C/f(T) + f/|(T))rfT. 



e«*yi(a;,t) = 5 + (e«^ - l)e + Ce^ + Ke^^Yx{T)dT + CJ^ e«^(5f (r) + Yi{T))dT. 



e''*Y2{x, t)= S+ {e^'^ - l)e + Ce^ + Ke''^y2(T)dT + C/; e'^^{Y^{T) + Yi{T))dT. 
Then obviously, Ui < Yi for i = 1. 2 and Yi satisfies 



(4.26) 

Therefore W te[T,T + Ti], 
(4.27) 



dt 



Yi{t)=C{Y,' + Y^'), 



r,(T) = e-''^ S + (e''^ - l)s + Ce 



Y{t) < e-'^'^' + {e"''' - l)s + Cs'^ + C J^Yirfdr, 
where Y{t) = {Yi{t),Y2{t))'^ . More precisely, we have V t e [T,T + Ti], 

(4.28) Y{t) < 





(5 + (e-*^ - l)e - 


f 




1 - Ce-«^ 


(5 + (e«^ - l)e4 


-C£2 


(t-T) 



Therefore for any given Tq > 0, we can choose < £o < 1/ (2C) sufficiently small and independent of T. For 
any fixed s {0 < e < eo), there exists 5 = 6{s) > also independent of T, such that 



(4.29) 



yit)< 



e-«^5 + (1 - e-'^^/2)£ 



1-C 



e-«T^+(l-e-«r/2)e 



<£, e [r,r + ri]. 



Therefore we may choose Eq so small that 2C£o < 1 and CTiEq < l/(4e'*"^'' — 1). Then choose S < e/A. In this 

way, we further have Tq ~ | log£|. This proves the theorem. 

□ 
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5. Conclusion. We studied the free boundary problem of swelling of a gel with inertia and elastic effects 

dominating the viscous ones, in one space dimension. This corresponds to early dynamics regime for many 
polymeric gels and also describes the main dynamics for polyssaccharide gel networks. We transform the free 
boundary problem of the weakly dissipative hyperbolic governing equation to one fixed domain. As a result of 
requiring the system to admit one dimensional solutions, it turns out that the polymer volume fraction at the 
interface with the surrounding fluid is fully determined, if we neglect fluid inertia in the boundary conditions. 

We find necessary and sufficient conditions for the system to be hyperbolic, this requiring a monotonicity 
property of a one-dimensional stress. In current work, we are studying the case that such condition fails, and 
found that it may be indicative of two possible phenomena, one being de-swelling, that is the interface will 
move backwards, and the second one may be related to the gel locally experiencing a volume phase transition. 
This is a phenomenon analogous to the volume transition in polyelectrolyte gels that occur when a relevant 
ion concentration reaches a critical value. 
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